Uncertainty visualization using HDR volume rendering

In this paper, we explore a novel idea of using high dynamic range (HDR) technology for uncertainty visualization. We focus on scalar volumetric data sets where every data point is associated with scalar uncertainty. We design a transfer function that maps each data point to a color in HDR space. The luminance component of the color is exploited to capture uncertainty. We modify existing tone mapping techniques and suitably integrate them with volume ray casting to obtain a low dynamic range (LDR) image. The resulting image is displayed on a conventional 8-bits-per-channel display device. The usage of HDR mapping reveals fine details in uncertainty distribution and enables the users to interactively study the data in the context of corresponding uncertainty information. We demonstrate the utility of our method and evaluate the results using data sets from ocean modeling.

Modeling Indian Ocean circulation: Bay of Bengal fresh water plume and the Arabian Sea Mini Warm Pool

The Indian subcontinent divides the north Indian Ocean into two tropical basins, namely the Arabian Sea and the Bay of Bengal. The Arabian Sea has high salinity whereas the salinity of the Bay of Bengal is much lower due to the contrast in freshwater forcing of the two basins. The freshwater received by the Bay in large amounts during the summer monsoon through river discharge is flushed out annually by ocean circulation. After the withdrawal of the summer monsoon, the Ganga – Brahmaputra river plume flows first along the Indian coast and then around Sri Lanka into the Arabian Sea creating a low salinity pool in the southeastern Arabian Sea (SEAS). In the same region, during the pre-monsoon months of February – April, a warm pool, known as the Arabian Sea Mini Warm Pool (ASMWP), which is distinctly warmer than the rest of the Indian Ocean, takes shape. In fact, this is the warmest region in the world oceans during this period. Simulation of the river plume and its movement as well as its implications to thermodynamics has been a challenging problem for models of Indian Ocean. Here we address these issues using an ocean general circulation model – first we show that the model is capable of reproducing fresh plumes in the Bay of Bengal as well as its movement and then we use the model to determine the processes that lead to formation of the ASMWP. Hydrographic observations from the western Bay of Bengal have shown the presence of a fresh plume along the northern part of the Indian coast during summer monsoon. The Indian Ocean model when forced by realistic winds and climatological river discharge reproduces the fresh plume with reasonable accuracy. The fresh plume does not advect along the Indian coast until the end of summer monsoon. The North Bay Monsoon Current, which flows eastward in the northern Bay, separates the low salinity water from the more saline southern parts of the bay and thus plays an important role in the fresh water budget of the Bay of Bengal. The model also reproduces the surge of the fresh-plume along the Indian coast, into the Arabian Sea during northeast monsoon. Mechanisms that lead to the formation of the Arabian Sea Mini Warm Pool are investigated using several numerical experiments. Contrary to the existing theories, we find that salinity effects are not necessary for the formation of the ASMWP. The orographic effects of the Sahyadris (Western Ghats) and resulting reduction in wind speed leads to the formation of the ASMWP. During November – April, the SEAS behave as a low-wind heatdominated regime where the evolution of sea surface temperature is solely determined by atmospheric forcing. In such regions the evolution of surface layer temperature is not dependent on the characteristics of the subsurface ocean such as the barrier layer and temperature inversion.

Longwave radiative forcing of Indian Ocean tropospheric aerosol

A spectrally resolved discrete-ordinates radiative transfer model is used to calculate the change in downwelling surface and top-of-the-atmosphere (TOA) outgoing longwave (3.9-500 mum) radiative fluxes induced by tropospheric aerosols of the type observed over the Indian Ocean during the Indian Ocean Experiment (INDOEX). Both external and internal aerosol mixtures were considered. Throughout the longwave, the aerosol volume extinction depends more strongly on relative humidity than in most of the shortwave (0.28-3.9 mum), implying that particle growth factors and realistic relative humidity profiles must be taken into account when modeling the longwave radiative effects of aerosols. A typical boundary layer aerosol loading, with a 500-nm optical depth of 0.3, will increase the downwelling longwave flux at the surface by 7.7 W m(-2) over the clean air case while decreasing the outgoing longwave radiation by 1.3 W m(-2). A more vertically extended aerosol loading, exhibiting a high opacity plume between 2 and 3 km above the surface and having a typical 500-nm optical depth of 0.7, will increase the downwelling longwave flux at the surface by 11.2 W m(-2) over the clean air case while decreasing the outgoing longwave radiation by 2.7 W m(-2). For a vertically extended aerosol profile, approximately 30% of the TOA radiative forcing comes from sea salt and approximately 60% of the forcing comes from the combination of sea salt and dust. The remaining forcing is from anthropogenic constituents. These results are for the external mixture. For an internal mixture, TOA longwave forcings can be up to a factor of two larger. Therefore, to complete our understanding of this region's longwave aerosol radiative properties, more detailed information is needed about aerosol mixing states. These longwave radiative effects partially offset the large shortwave aerosol radiative forcing and should be included in regional and global climate modeling simulations.

Modeling the vulnerability of an urban groundwater system due to the combined impacts of climate change and management scenarios

Climate change impact on a groundwater-dependent small urban town has been investigated in the semiarid hard rock aquifer in southern India. A distributed groundwater model was used to simulate the groundwater levels in the study region for the projected future rainfall (2012-32) obtained from a general circulation model (GCM) to estimate the impacts of climate change and management practices on groundwater system. Management practices were based on the human-induced changes on the urban infrastructure such as reduced recharge from the lakes, reduced recharge from water and wastewater utility due to an operational and functioning underground drainage system, and additional water extracted by the water utility for domestic purposes. An assessment of impacts on the groundwater levels was carried out by calibrating a groundwater model using comprehensive data gathered during the period 2008-11 and then simulating the future groundwater level changes using rainfall from six GCMs Institute of Numerical Mathematics Coupled Model, version 3.0 (INM-CM. 3.0); L'Institut Pierre-Simon Laplace Coupled Model, version 4 (IPSL-CM4); Model for Interdisciplinary Research on Climate, version 3.2 (MIROC3.2); ECHAM and the global Hamburg Ocean Primitive Equation (ECHO-G); Hadley Centre Coupled Model, version 3 (HadCM3); and Hadley Centre Global Environment Model, version 1 (HadGEM1)] that were found to show good correlation to the historical rainfall in the study area. The model results for the present condition indicate that the annual average discharge (sum of pumping and natural groundwater outflow) was marginally or moderately higher at various locations than the recharge and further the recharge is aided from the recharge from the lakes. Model simulations showed that groundwater levels were vulnerable to the GCM rainfall and a scenario of moderate reduction in recharge from lakes. Hence, it is important to sustain the induced recharge from lakes by ensuring that sufficient runoff water flows to these lakes.

Copula-Based Modeling of TMI Brightness Temperature With Rainfall Type

Overland rain retrieval using spaceborne microwave radiometer offers a myriad of complications as land presents itself as a radiometrically warm and highly variable background. Hence, land rainfall algorithms of the Tropical Rainfall Measuring Mission (TRMM) Microwave Imager (TMI) have traditionally incorporated empirical relations of microwave brightness temperature (Tb) with rain rate, rather than relying on physically based radiative transfer modeling of rainfall (as implemented in the TMI ocean algorithm). In this paper, sensitivity analysis is conducted using the Spearman rank correlation coefficient as benchmark, to estimate the best combination of TMI low-frequency channels that are highly sensitive to the near surface rainfall rate from the TRMM Precipitation Radar (PR). Results indicate that the TMI channel combinations not only contain information about rainfall wherein liquid water drops are the dominant hydrometeors but also aid in surface noise reduction over a predominantly vegetative land surface background. Furthermore, the variations of rainfall signature in these channel combinations are not understood properly due to their inherent uncertainties and highly nonlinear relationship with rainfall. Copula theory is a powerful tool to characterize the dependence between complex hydrological variables as well as aid in uncertainty modeling by ensemble generation. Hence, this paper proposes a regional model using Archimedean copulas, to study the dependence of TMI channel combinations with respect to precipitation, over the land regions of Mahanadi basin, India, using version 7 orbital data from the passive and active sensors on board TRMM, namely, TMI and PR. Studies conducted for different rainfall regimes over the study area show the suitability of Clayton and Gumbel copulas for modeling convective and stratiform rainfall types for the majority of the intraseasonal months. Furthermore, large ensembles of TMI Tb (from the most sensitive TMI channel combination) were generated conditional on various quantiles (25th, 50th, 75th, and 95th) of the convective and the stratiform rainfall. Comparatively greater ambiguity was observed to model extreme values of the convective rain type. Finally, the efficiency of the proposed model was tested by comparing the results with traditionally employed linear and quadratic models. Results reveal the superior performance of the proposed copula-based technique.

Slow Slip Acceleration beneath Andaman Islands Triggered by the 11 April 2012 Indian Ocean Earthquakes

The M-w 8.6 and 8.2 strike-slip earthquakes that struck the northeast Indian Ocean on 11 April 2012 resulted in coseismic deformation both at near and distant sites. The slip distribution, deduced using seismic-wave analysis for the orthogonal faults that ruptured during these earthquakes, is sufficient to predict the coseismic displacements at the Global Positioning System (GPS) sites, such as NTUS, PALK, and CUSV, but fall short at four continuous sites in the Andaman Islands region. Slip modeling, for times prior to the events, suggests that the lower portion of the thrust fault beneath the Andaman Islands has been slipping at least at the rate of 40 cm/yr, in response to the 2004 Sumatra-Andaman coseismic stress change. Modeling of GPS displacements suggests that the en echelon and orthogonal fault ruptures of the 2012 intraplate oceanic earthquakes could have possibly accelerated the ongoing slow slip, along the lower portion of the thrust fault beneath the islands with a month-long slip of 4-10 cm. The misfit to the coseismic GPS displacements along the Andaman Islands could be improved with a better source model, assuming that no local process contributed to this anomaly.

A modeling study of the processes of surface salinity seasonal cycle in the Bay of Bengal

In response to the Indian Monsoon freshwater forcing, the Bay of Bengal exhibits a very strong seasonal cycle in sea surface salinity (SSS), especially near the mouths of the Ganges-Brahmaputra and along the east coast of India. In this paper, we use an eddy-permitting (similar to 25 km resolution) regional ocean general circulation model simulation to quantify the processes responsible for this SSS seasonal cycle. Despite the absence of relaxation toward observations, the model reproduces the main features of the observed SSS seasonal cycle, with freshest water in the northeastern Bay, particularly during and after the monsoon. The model also displays an intense and shallow freshening signal in a narrow (similar to 100 km wide) strip that hugs the east coast of India, from September to January, in good agreement with high-resolution measurements along two ships of opportunity lines. The mixed layer salt budget confirms that the strong freshening in the northern Bay during the monsoon results from the Ganges-Brahmaputra river discharge and from precipitation over the ocean. From September onward, the East India Coastal Current transports this freshwater southward along the east coast of India, reaching the southern tip of India in November. The surface freshening results in an enhanced vertical salinity gradient that increases salinity of the surface layer by vertical processes. Our results reveal that the erosion of the freshwater tongue along the east coast of India is not driven by northward horizontal advection, but by vertical processes that eventually overcome the freshening by southward advection and restore SSS to its premonsoon values. The salinity-stratified barrier layer hence only acts as a ``barrier'' for vertical heat fluxes, but is associated with intense vertical salt fluxes in the Bay of Bengal.

Joint modeling of lithosphere and mantle dynamics: Evaluation of constraints from global tomography models

With the advances in technology, seismological theory, and data acquisition, a number of high-resolution seismic tomography models have been published. However, discrepancies between tomography models often arise from different theoretical treatments of seismic wave propagation, different inversion strategies, and different data sets. Using a fixed velocity-to-density scaling and a fixed radial viscosity profile, we compute global mantle flow models associated with the different tomography models and test the impact of these for explaining surface geophysical observations (geoid, dynamic topography, stress, and strain rates). We use the joint modeling of lithosphere and mantle dynamics approach of Ghosh and Holt (2012) to compute the full lithosphere stresses, except that we use HC for the mantle circulation model, which accounts for the primary flow-coupling features associated with density-driven mantle flow. Our results show that the seismic tomography models of S40RTS and SAW642AN provide a better match with surface observables on a global scale than other models tested. Both of these tomography models have important similarities, including upwellings located in Pacific, Eastern Africa, Iceland, and mid-ocean ridges in the Atlantic and Indian Ocean and downwelling flows mainly located beneath the Andes, the Middle East, and central and Southeast Asia.

Modeling of Channel Potential and Subthreshold Slope of Symmetric Double-Gate Transistor

A new physically based classical continuous potential distribution model, particularly considering the channel center, is proposed for a short-channel undoped body symmetrical double-gate transistor. It involves a novel technique for solving the 2-D nonlinear Poisson's equation in a rectangular coordinate system, which makes the model valid from weak to strong inversion regimes and from the channel center to the surface. We demonstrated, using the proposed model, that the channel potential versus gate voltage characteristics for the devices having equal channel lengths but different thicknesses pass through a single common point (termed ``crossover point''). Based on the potential model, a new compact model for the subthreshold swing is formulated. It is shown that for the devices having very high short-channel effects (SCE), the effective subthreshold slope factor is mainly dictated by the potential close to the channel center rather than the surface. SCEs and drain-induced barrier lowering are also assessed using the proposed model and validated against a professional numerical device simulator.

Modeling of sonochemical decomposition of CCI4 in aqueous solutions

In the context of removal of organic pollutants from wastewater, sonolysis of CCl4 dissolved in water has been widely investigated. These investigations are either completely experimental or correlate data empirically. In this work, a quantitative model is developed to predict the rate of sonolysis of aqueous CCl4. The model considers the isothermal growth and partially adiabatic collapse of cavitation bubbles containing gas and vapor leading to conditions of high temperatures and pressures in them, attainment of thermodynamic equilibrium at the end of collapse, release of bubble contents into the liquid pool, and reactions in the well-mixed pool. The model successfully predicts the extent of degradation of dissolved CCl4, and the influence of various parameters such as initial concentration of CCl4, temperature, and nature of gas atmosphere above the liquid. in particular, it predicts the results of Hua and Hoffmann (Environ. Sci Technol, 1996, 30, 864-871), who found that degradation is first order with CCl4 and that Argon as well as Ar-O-3 atmospheres give the same results. The framework of the model is capable of quantitatively describing the degradation of many dissolved organics by considering all the involved species.

Multiscale modeling of lamellar mesophases

The mesoscale simulation of a lamellar mesophase based on a free energy functional is examined with the objective of determining the relationship between the parameters in the model and molecular parameters. Attention is restricted to a symmetric lamellar phase with equal volumes of hydrophilic and hydrophobic components. Apart from the lamellar spacing, there are two parameters in the free energy functional. One of the parameters, r, determines the sharpness of the interface, and it is shown how this parameter can be obtained from the interface profile in a molecular simulation. The other parameter, A, provides an energy scale. Analytical expressions are derived to relate these parameters to r and A to the bending and compression moduli and the permeation constant in the macroscopic equation to the Onsager coefficient in the concentration diffusion equation. The linear hydrodynamic response predicted by the theory is verified by carrying out a mesoscale simulation using the lattice-Boltzmann technique and verifying that the analytical predictions are in agreement with simulation results. A macroscale model based on the layer thickness field and the layer normal field is proposed, and the relationship between the parameters in the macroscale model from the parameters in the mesoscale free energy functional is obtained.

Molecule Modeling of Human Pentameric alpha(7) Neuronal Nicotinic Acetylcholine Receptor and Its Interaction with its Agonist and Competitive Antagonist

The nicotinic Acetylcholine Receptor (nAChR) is the major class of neurotransmitter receptors that is involved in many neurodegenerative conditions such as schizophrenia, Alzheimer's and Parkinson's diseases. The N-terminal region or Ligand Binding Domain (LBD) of nAChR is located at pre- and post-synaptic nervous system, which mediates synaptic transmission. nAChR acts as the drug target for agonist and competitive antagonist molecules that modulate signal transmission at the nerve terminals. Based on Acetylcholine Binding Protein (AChBP) from Lymnea stagnalis as the structural template, the homology modeling approach was carried out to build three dimensional model of the N-terminal region of human alpha(7)nAChR. This theoretical model is an assembly of five alpha(7) subunits with 5 fold axis symmetry, constituting a channel, with the binding picket present at the interface region of the subunits. alpha-netlrotoxin is a potent nAChR competitive antagonist that readily blocks the channel resulting in paralysis. The molecular interaction of alpha-Bungarotoxin, a long chain alpha-neurotoxin from (Bungarus multicinctus) and human alpha(7)nAChR seas studied. Agonists such as acetylcholine, nicotine, which are used in it diverse array of biological activities, such as enhancements of cognitive performances, were also docked with the theoretical model of human alpha(7)nAChR. These docked complexes were analyzed further for identifying the crucial residues involved i interaction. These results provide the details of interaction of agonists and competitive antagonists with three dimensional model of the N-terminal region of human alpha(7)nAChR and thereby point to the design of novel lead compounds.

Analysis Of The Energy Quantization Effects On Single Electron Inverter Performance Through Noise Margin Modeling

Possible integration of Single Electron Transistor (SET) with CMOS technology is making the study of semiconductor SET more important than the metallic SET and consequently, the study of energy quantization effects on semiconductor SET devices and circuits is gaining significance. In this paper, for the first time, the effects of energy quantization on SET inverter performance are examined through analytical modeling and Monte Carlo simulations. It is observed that the primary effect of energy quantization is to change the Coulomb Blockade region and drain current of SET devices and as a result affects the noise margin, power dissipation, and the propagation delay of SET inverter. A new model for the noise margin of SET inverter is proposed which includes the energy quantization effects. Using the noise margin as a metric, the robustness of SET inverter is studied against the effects of energy quantization. It is shown that SET inverter designed with CT : CG = 1/3 (where CT and CG are tunnel junction and gate capacitances respectively) offers maximum robustness against energy quantization.

Ocean atmosphere coupling over monsoon regions

In monsoon regions, the seasonal migration of the intertropical convergence zone (ITCZ) is manifested as a seasonal reversal of winds. Most of the summer monsoon rainfall over India occurs owing to synoptic and large-scale convection associated with the continental ITCZ (Fig. 1). We have investigated the interaction between these large-scale convective systems and the ocean over which they are generated1â3, concentrating on the relationship between organized convection over the Indian Ocean and sea surface temperature (SST). We report here that on a monthly basis the degree of cloudiness correlates well with SST for the relatively colder oceans, but when SST is maintained above 28 °C it ceases to be an important factor in determining the variability of cloudiness. Over the major regions of convection east of 70°E, which are warm year after year, the observed cloudiness cannot be correlated with variations in SST.

Modeling of a pressure modulated desalination system using bond graph methodology

A desalination system is a complex multi energy domain system comprising power/energy flow across several domains such as electrical, thermal, and hydraulic. The dynamic modeling of a desalination system that comprehensively addresses all these multi energy domains is not adequately addressed in the literature. This paper proposes to address the issue of modeling the various energy domains for the case of a single stage flash evaporation desalination system. This paper presents a detailed bond graph modeling of a desalination unit with seamless integration of the power flow across electrical, thermal, and hydraulic domains. The paper further proposes a performance index function that leads to the tracking of the optimal chamber pressure giving the optimal flow rate for a given unit of energy expended. The model has been validated in steady state conditions by simulation and experimentation.